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Abstract 

This  paper  considers  a  connection  between  the  deterministic  and  noisy  behavior 
of  nonlinear  networks.  Specifically,  a  particular  bridge  circuit  is  examined  which  has 
two  possibly  nonlinear  energy  storage  elements.  By  proper  choice  of  the  constitutive 
relations  for  the  network  elements,  the  deterministic  terminal  behavior  reduces  to  that 
of  a  single  linear  resistor.  This  reduction  of  the  deterministic  terminal  behavior,  in 
which  a  natural  frequency  of  a  linear  circuit  does  not  appear  in  the  driving-point 
impedance,  has  been  shown  in  classical  circuit  theory  books  {e.g.  mm)-  The  paper 
shows  that,  in  addition  to  the  reduction  of  the  deterministic  behavior,  the  thermal 
noise  at  the  terminals  of  the  network,  arising  from  the  usual  Nyquist- Johnson  noise 
model  associated  with  each  resistor  in  the  network,  is  also  exactly  that  of  a  single 
linear  resistor.  While  this  result  for  the  linear  time-invariant  (LTI)  case  is  a  direct 
consequence  of  a  well-known  result  for  RLC  circuits,  the  nonlinear  result  is  novel.  We 
show  that  the  terminal  noise  current  is  precisely  that  predicted  by  the  Nyquist-Johnson 
model  for  R  if  the  driving  voltage  is  zero  or  constant,  but  not  if  the  driving  voltage  is 
time-dependent  or  the  inductor  and  capacitor  are  time-varying. 
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Figure  1:  The  linear  noise-free  bridge  circuit  is  matched  and  has  input  impedance  R  is 
L  =  R^C. 

1  Introduction 

Consider  the  bridge  circuit  of  Figure  [Tj  It  is  a  standard  result  of  linear  circuit  theory  that 
under  the  matching  condition  L  =  R^  C  the  natural  frequency  of  the  circuit  does  not  appear 
as  a  pole  in  the  driving-point  impedance  mi.  One  can  intuitively  see  that  -  regardless  of 
the  values  of  the  capacitor  and  inductor  -  for  high  frequencies,  the  capacitor  is  essentially 
a  short  circuit,  whereas  the  inductor  is  essentially  an  open  circuit;  at  low  frequencies,  the 
opposite  occurs.  The  matching  condition  ensures  that  a  balance  is  preserved  for  intermediate 
frequencies:  the  charging  of  the  capacitor  is  matched  by  the  fluxing  of  the  inductor. 

Example  2  on  pp.  630-633  of  Ref.  pQ  discusses  the  circuit  shown  in  Figure]^  This  circuit 
will  be  the  main  focus  of  this  paper.  However,  the  derivations  can  also  be  carried  out  for  the 
dual  circuit,  which  is  in  fact  the  version  presented  in  Ref.  |2].  Ref.  |2]  also  discusses  the  circuit 
in  terms  of  loss  of  observability  and  controllability:  under  the  matching  condition,  the  states 
of  the  system,  namely  the  inductor  flux  and  the  capacitor  charge,  are  neither  observable  nor 
controllable  from  the  external  terminals,  i.e.,  the  state  equations  become  nonminimal. 

It  is  straightforward  to  verify  directly  in  this  linear  case  that  if  a  Nyquist- Johnson  Elll] 
noise  model  (as  seen  in  Fig.  is  associated  with  each  resistor,  then  the  spectrum  of  the 
short-circuit  current  is  also  that  predicted  by  a  Nyquist- Johnson  noise  model  for  a  single 
resistor  of  value  R.  The  verihcation  can  be  done  by  standard  frequency-domain  techniques 
or  by  stochastic  calculus.  The  highpass  hltering  of  the  RC  branch  is  precisely  balanced 
by  the  lowpass  hltering  of  the  RL  branch,  so  that  the  terminal  noise  spectrum  is  hat.  Of 
course,  both  resistors  must  be  at  the  same  temperature.  As  noted  in  0,  applying  a  d.c. 
voltage  to  the  circuit  would  result  in  diherential  heating  of  the  resistor  in  the  RL  branch. 
If  the  resistors  were  not  properly  connected  to  thermal  reservoirs,  one  could  heat  up  and 
become  noisier  than  the  other,  and  the  noise  spectrum  would  no  longer  be  hat.  This  is  a 
trivial  nonequilibrium  exception  to  the  results  of  this  paper,  which  assumes  uniform,  constant 
temperature. 

The  result  above  is  a  particular  example  of  a  general  circuit  theory  result,  namely,  that  a 
one-port  network  of  (linear)  passive  elements  with  port  impedance  Y [jut)  presents  a  thermal 
noise  voltage  with  power  spectrum  2kT  Re{y(ja;)},  where  k  is  Boltzmann’s  constant  and 
T  is  the  absolute  temperature  [3].  Physicists  regard  such  results  as  particular  cases  of  the 
huctuation-dissipation  theorem  [6]. 
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Figure  2:  Nyquist-Johnson  thermal  noise  model  (Norton  form)  is  a  noiseless  linear  resistor  in 
parallel  with  a  Gaussian  white  noise  current  source  with  power  spectral  density  2kT/R. 


This  paper  extends  these  results  to  the  nonlinear  case  and  to  some  nonequilibrium  situ¬ 
ations.  In  layman’s  terms,  the  question  is  as  follows:  suppose  one  has  two  black  boxes,  one 
with  a  matched  bridge  inside,  and  one  with  a  single  equivalent  linear  resistor.  Is  it  possible 
to  distinguish  the  two  using  the  noise  behavior? 

This  paper  studies  one  carefully-chosen  example,  motivated  by  the  question  of  whether 
some  form  of  fluctuation-dissipation  theorem  holds  for  some  class  of  nonlinear  circuits.  Our 
initial  formulation  appears  below  as  a  conjecture  for  any  pair  of  two-terminal  networks,  each 
comprising  an  interconnection  of  LTl  resistors  at  a  uniform,  constant  temperature,  described 
by  the  Nyquist-Johnson  model,  and  possibly  also  capacitors  and  inductors  that  may  be 
nonlinear  or  time-varying.  Two  such  networks  are  said  to  be  zero-state  deterministically 
equivalent  if  every  applied  terminal  voltage  waveform  v{t),t  >  0,  produces  the  same  current 
response  i{t)  from  both  networks,  provided  all  capacitor  voltages  and  inductor  currents  are 
initially  zero  and  all  noise  sources  in  the  resistor  models  are  set  to  zero.  (In  the  LTl  case 
this  just  means  the  two  input  admittances  are  identical.) 

Preliminary  Fluctuation- Dissipation  Conjecture  for  Networks: 

No  two  zero-state  deterministically  equivalent  networks  can  be  distinguished  by 
their  terminal  noise  current  responses  to  any  applied  voltage  waveform. 

The  conjecture  just  hypothesizes  that  the  deterministic  terminal  behavior  uniquely  de¬ 
termines  the  noise  current  response  for  all  voltage  drives,  independent  of  the  details  of  the 
network.  The  conjecture  is  true  in  the  LTl  case.  (Closely  related  formulations  for  the  current- 
driven  and  multiport  cases  [5]  also  hold  true  for  LTl  networks,  but  we  ignore  them  here  for 
simplicity.)  However,  an  examination  of  the  bridge  circuit  will  show  that  this  conjecture  is 
wrong  in  other  cases  and  must  be  revised. 

The  linear  results  are  hrst  derived  in  Sections  2-5  to  motivate  the  nonlinear  case  and 
to  familiarize  the  reader  with  the  notation.  In  Section  6  we  derive  the  matching  condition 
corresponding  to  L  =  C  for  the  bridge  circuit  with  nonlinear,  time-invariant  inductor 
and  capacitor,  under  which  it  becomes  deterministically  equivalent  at  the  terminals  to  a 
single  linear  resistor  R.  In  Section  7,  we  show  that  such  a  matched  nonlinear  bridge  gives  a 
short-circuit  port  current  noise  statistically  identical  to  that  of  the  Nyquist-Johnson  model 
for  R  at  thermal  equilibrium.  Nonequilibrium  situations  are  considered  in  Section  8. 
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2  Linear,  Noise-Free  Case 


Kirchoff’s  Laws  for  the  bridge  circuit  of  Figure  [^yield  differential  equations  for  the  capacitor 
charge 


dq  _V  1 
dt~  R~ 


(1) 


and  for  the  inductor  current 


V-Ri. 


(2) 


Great  insight  can  be  extracted  from  the  quantity  i  —  q/RC.  Using  the  two  previous  equations, 
the  time  evolution  of  this  quantity  is 


dt  V  RC) 


RC^R  RC^ 


(3) 


Suppose  now  that  the  element  values  are  balanced  with  the  matching  condition 


Then,  Eq.  (|^  reduces  to 


L 

Wc 


1  or  L  =  R^C. 


dt  V  Rc)  L  V  Rc)  ■ 


(4) 

(5) 


With  zero  initial  conditions,  or  for  non-zero  initial  conditions  satisfying  i(0)  =  q{0)/RC,  this 
guarantees  that  i{t)  =  q{t)/RC  for  all  time.  More  signihcantly,  all  other  initial  conditions 
are  exponentially  attracted  to  the  line  i  =  q/RC.  Further,  by  Kirchoff’s  Current  Law, 


I 


dq  .  V  q 

V 

R' 


(6) 


so  that  the  entire  circuit  appears  as  simply  a  single  resistor  of  value  R.  The  mode  of  the 
circuit  corresponding  to  i  —  q/RC  is  neither  observable  (it  does  not  appear  in  the  output 
equation  for  I)  nor  controllable  (its  dynamics  are  unaffected  by  the  input  V). 
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I 


Figure  3:  Noisy  bridge  circuit,  Norton  form 


3  Linear  Noisy  Case 


To  analyze  the  noise  behavior  of  the  linear  bridge  circuit,  a  Norton  equivalent  current  noise 
source  is  associated  with  each  resistor,  as  in  Figure]^  The  noise  sources  ijs^  and  are  in¬ 
dependent,  white,  Gaussian  random  processes  with  spectral  density  2kT / R.  The  differential 
equations  for  the  circuit  become 


The  port  current  is 


dq  _  V  1 

~di  ~ 

di 

L-  =  V-R{i  +  iN,). 


^  dq  1/  1 


q-iN2- 


(7) 

(8) 

(9) 


The  goal  of  this  section  is  to  show  that  the  short-circuit  {V  =  0)  power  spectral  density 
for  I  is  2kT/R,  under  the  matching  condition  L/C  =  R^.  In  anticipation  of  the  nonlinear 
case,  where  frequency-domain  techniques  fail,  time-domain  techniques  are  used  to  hnd  the 
autocorrelation  function  for  I ;  the  Fourier  transform  then  yields  the  power  spectral  density. 
Since  there  are  no  explicit  time  dependencies,  the  random  process  I  is  stationary,  and  the 
autocorrelation  only  depends  on  the  difference  between  the  two  time  points. 


Rjj{t)  =  E{I{t)I{t  +  T)} 


=  E 

=  E 


Kt)  -  (lit)  - 


tit  + 


i(t)  i{t  +  r 
1 


RC 


i{t)  q{t  +  r)  -  i{t)  +  r) 


qit)  iit  +  r)  +  q{t)  q{t  +  ^)  +  ^  ^N^it  +  r) 

-iN2it)  ^(^  +  ^)  +  ;^  +  ^)  +  tN2it)  iN2it  +  t) 
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For  =  0,  the  two  branches  of  the  circuit  do  not  interact.  Since  the  two  noise  sources  are 
independent,  i(-)  is  not  correlated  with  either  g(-)  or  iN2{')-  For  t  >  0,  q(t)  is  uncorrelated 
with  +  t)  (the  capacitor  charge  now  is  not  affected  by  the  current  noise  in  the  future). 
Correspondingly,  for  r  <  0,  q(t  +  r)  is  uncorrelated  with  iN2(t)-  Since  each  of  the  quanti¬ 
ties  i(-),  g(-),  and  iN2{')  are  individually  zero-mean,  hve  of  the  nine  terms  in  the  equation 
immediately  vanish.  For  r  >  0, 

Rii{t)  =  E  +  +  ^^q{t)  q{t  +  T) +iN2{t)iN2{t  +  T) 

+  g(t  +  r)}  .  (10) 


The  time  arguments  of  the  last  term  will  be  switched  for  r  <  0. 

The  terms  in  this  expression  will  be  computed  individually.  First,  of  course, 

2kT 

E{^M2{t)^N2{t  +  r)}  =  —S{r). 


(11) 


The  two  branches  are  disjoint  for  V  =  0,  and  the  differential  equations  for  the  random 
variables  i  and  q  are  separable: 

di 

L—  =  -Ri-  RiMi, 


which  has  the  solution 


and 


which  has  the  solution 


i{t)  =  -  y  f  e  iNiit')  dt\ 

^  J  —OO 

dq  1 

J  —OO 


(12) 


(13) 


These  two  time-domain,  sample-path  solutions  are  the  key  to  computing  the  remaining 
expectations.  Assuming  zero  initial  conditions  (i(— cx))  =  0  and  q{—oo)  =  0)  so  that  the 
random  processes  are  stationary. 


E  {i{t)  i{t  +  r)}  =  E 


J 


R 


ft+T 


iNiit'')  dt'‘ 


td2  pt  pt+T 

=  —  /  /  e-d+^-^''WL  E{ij^^(t')iN,{f)}dt’ df 


=  —  C  ^-{t-t')RIL 


’  —OO 

pt+r 


2kT 


L2 

R^  2kT 


e-d+r-t")RIL 


—  OO 

'2  pt 


R 


^—{t—t'+t+T—t')R/L 
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2kTR 


^-i2t+r)R/L  /  ^2t'RIL 


,2t+T)R/L^  \p2tR/L  _  pi 

L2  2R^  J 


=  — 


In  fact,  the  above  derivation  assumed  r  >  0;  the  correct  expression  for  any  r  is 


Lj 


Similarly, 


( r  pt  If  pt+'T 

E  {q(t)  q(t  +  t)}  =  E  <  —  iN2(t')  dt'  —  dt” 

L  L  'J  —oo  J  L  'J  —oo 


^  ^-(2t+T)IRC 


rt  rt+T 


+t  E  iN2id'')}  dt'  dt" 


2kT  RC 
“R 


-\r\/RC 


=  kTCe-^^dRC_ 

The  last  term  to  be  computed  is 


E  {iN2{t)  q{t +  t)}  =  E  UN2{t)  -  j  e  iN2id')  dt' 

L  L  ^ — oo 


d^^E{iN2it)  iN2it')}dt' 


=  _^-(t+r)R/L  ^t'/RC 


=  _^-(t+T)IRC  ^t/RC 


=  r>0.  (17) 

R 

Unlike  the  previous  terms,  this  last  expectation  is  zero  for  r  <  0.  For  r  <  0,  the  time 
arguments  on  the  left-hand  side  switch  places,  resulting  in 

2kT 

E  +  t)  q{t)}  = - 
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Assembling  Eqs.  ([TT|,  (|l^,  (16),  and  ([TT])  and  substituting  into  Eq.  (10), 

1 


=  E  {i{t)  i{t  +  t)]  + 


E  {q{t)  q{t  +  t)}  +  E  (t)  In^ {t  +  t)} 


RC 


EiiN^it)  qit  +  r)} 


1  2kT 


-\t\/RC 


RC  R 


(19) 


Eq.  (18)  ensures  that  this  expression  is  valid  for  all  r. 

Under  the  matching  condition  L  =  R^C  among  the  element  values,  the  leading  coefficients 
are  related: 


kT  _  kT 

~T  ~  Wc' 

and  the  time  constants  in  the  exponential  are  equal: 

1  _  R 

Wj~T 


Therefore,  Eq.  (19)  simplihes  to 


9kT 

Rnir)  =  ^(r). 


(20) 


(21) 


(22) 


The  Fourier  transform  of  a  delta-function  is  simple  to  compute,  immediately  yielding  the 
desired  result, 

2kT 

Sjj{u;)=E{RniT)}  =  —,  (23) 

for  all  frequencies  oj.  The  noise  behavior  of  the  circuit  is  equivalent  to  that  of  a  single  linear 
resistor,  just  as  was  the  case  for  the  deterministic  behavior. 

This  same  result  may  also  be  obtained  without  short-circuiting  the  terminals.  If  the 
terminals  are  driven  by  a  deterministic  voltage  waveform  V{t),  the  response  of  the  matched 
bridge  circuit  will  be  indistinguishable  from  the  response  of  a  single  linear  resistor.  The 
autocorrelation  of  the  output  current  will  reduce  to 


R„{t)  =  E  {/(()  I{t  +  t)}=  i(T), 


R 


(24) 


the  deterministic  autocorrelation  of  the  applied  voltage  driving  current  through  a  linear 
resistor,  plus  the  autocorrelation  of  the  noise  of  a  single  linear  resistor.  This  fact  will  be 
proved  in  Appendix  I. 


Figure  4:  Linear  bridge  circuit  in  the  frequency  domain 


4  Frequency  Domain 

The  result  of  Section  can  also  be  obtained  using  standard  frequency-domain  techniques  for 
the  power  spectral  density.  The  Norton-form  noise  sources  drive  an  RC  or  an  RL  hlter  and 
the  responses  are  combined  to  form  the  output  noise  current. 

First,  the  matching  condition  will  be  re-derived  using  the  complex  impedance  in  the 
frequency  domain.  Let  be  the  impedance  of  the  RL  branch  of  Figure 

Z\  =  Z^  -|-  Zi^  =  R  +  jujL,  (25) 

where  j  =  \/—l,  and  let  Z2  be  the  impedance  of  the  RC  branch: 

Z2  =  Zr  +  Zc  =  R  +  1/ (jcvC).  (26) 


The  port  impedance  Z  =  V/ 1  is  given  by  the  parallel  combination  of  these  two  impedances 

^  ^  Z1Z2  {R  +  ju;L){R+l/{ju;C)) 

Zi  -\-  Z2  R  jujL  -1-  i?  -1- 1/ 

R^  +  R/  {jujC)  +JUJLR  + L/C 
2R  +  juL  -1-1/ (jujC) 

R  +  L/iRC)+ju;L  +  l/{ju;C) 

2R  +  ju;L  +  l/{juC)  '  ^  ^ 

Clearly,  the  fraction  is  unity  if  L/{RC)  =  R,  which  is  the  matching  condition. 

Now,  associate  with  each  resistor  a  Norton-form  Gaussian  white-noise  current  source 
with  power  spectral  density  2kT / R.  The  sources  will  be  connected  to  input  ports  1  and  2  in 
Figure]^  Because  of  the  short-circuit  condition,  the  bridge  can  be  broken  into  two  parts,  as 
shown,  for  the  purposes  of  computing  the  current  due  to  each  source.  The  output  current  is 
a  combination  of  these  two  currents,  and  since  the  two  noise  sources  are  independent,  their 
respective  contributions  to  the  output  current  power  spectral  density  add. 
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For  the  noise  source  connected  to  port  1,  its  current  is  divided  between  the  resistor  and 
inductor,  and  its  contribution  to  the  output  current  is  only  the  current  passing  through  the 
inductor.  If  a  hlter  H{uj)  is  applied  to  an  input  X{uj),  the  output  will  have  power  spectral 
density  Y{u)  =  \H{u)\‘^X{u).  The  hlter  in  this  case  has  the  shape  Zji/{Zji  +  Z^)  as  dictated 
by  current  division.  Thus,  the  output  current  power  spectral  density  for  Ji  is 


2kT 

Zr 

^  2kT 

R 

^  2kT  R^  2kTR 

R 

Zr  +  Zl 

R 

R  +  jujL 

~  R  R^  +  uj^L‘^  “  +  o;2L2 

(28) 


For  the  noise  source  connected  to  port  2,  similarly. 


2kT 

Zr 

^  2kT 

R 

^  2kTR 

R 

Zr  +  Zc 

R 

R  +  1/{3uC) 

“  +  l/(o;2C'2) 

(29) 


Since  the  two  noise  sources  are  independent,  the  output  voltage  power  spectral  density  is 
simply  the  sum  of  the  two  half-circuit  power  spectral  densities: 


Sii{uj) 


2kTR  ^  2kTR 


2kTR 


2kTR 


R^  +  +  R^  +  uj^L^\ 

[R^  +  u^L^][R^  +  l/{u^C^)]  ) 

2R^  +  l/{u^C^)+u^L^ 

+  Ry(uj^C^)  +  R^oj^L^  +  L^jC^ 


2kTR  /  +  l/(o;2C'2)  + 

R^  +  L7(i?2C'2)  +  l/(o^2C'2)  +  ^2^2 


2kT 


(30) 


where  the  last  equality  follows  because  the  fraction  inside  the  large  parentheses  on  the 
previous  line  is  unity  under  the  matching  condition  L‘^/{R^C‘^)  =  R^ . 
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5  Dual  Circuit 

The  circuit  of  Figure  is  the  dual  of  the  circuit  in  Fig.  [T]  the  parallel  combination  of 
two  branches,  each  containing  two  elements  in  series,  has  become  a  series  combination  of 
two  loops,  each  containing  two  elements  in  parallel.  Resistors  are  their  own  duals,  and  the 
inductor  and  capacitor  are  each  other’s  dual.  Lastly,  for  analyzing  the  noise,  the  open- 
circuit  voltage  noise  takes  the  place  of  the  short-circuit  current  noise.  The  construction  of  a 
dual  circuit  is  best  explained  in  Ref.  [1],  chapter  10,  section  4.  One  can  also  construct  the 
dual  circuit  for  the  nonlinear  bridge  presented  later  in  this  paper;  for  the  duals  of  nonlinear 
elements,  the  reader  is  referred  to  Ref. 

The  differential  equation  for  the  capacitor  charge  is 


dq  _  j  q 

dt  RC' 


(31) 


and  for  the  inductor  flux 

(The  state  variable  for  the  inductor  has  switched  to  flux  in  this  section  in  preparation  for 
the  nonlinear  sections  to  follow,  in  which  flux  is  a  more  natural  state  variable.)  These  follow 
directly  from  the  constitutive  relations  of  the  elements  and  Kirchoff’s  Laws.  Consider  the 
quantity  0/L  —  q/ RC.  Using  the  two  previous  equations,  the  time  evolution  of  this  quantity 

is 


dt 


Using  the  same  matching  condition  as  before,  namely,  ^ 


1 


Eq.  (33)  reduces  to 


jJ] 

dt  \L  RC J  L\L  RC J  ' 


(33) 


(34) 
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Corresponding  to  the  line  i  =  q/RC  in  the  original  circuit,  trajectories  for  the  dual  circuit 
are  exponentially  attracted  to  the  line  0/L  =  q/RC  in  the  (0,  g)  plane,  and  of  course  for 
any  initial  conditions  on  the  line  (including  the  origin),  the  trajectory  stays  on  the  line  for 
all  time.  Further,  by  Kirchoff’s  Voltage  Law, 


V 


/L 

C  dt 


1 

c 


=  IR, 


(35) 


so  that  the  entire  dual  circuit  also  appears  as  simply  a  single  resistor  of  value  R. 

The  dual  of  the  short-circuit  output  current  is  the  open-circuit  voltage.  The  open-circuit 
voltage  power  spectral  density  could  be  computed  by  either  of  the  previous  two  methods, 
namely,  via  the  expectation  over  sample  paths  or  by  frequency-domain  hltering.  Instead, 
a  third  technique  will  be  used  in  preparation  for  the  nonlinear  cases  of  later  chapters.  In 
the  nonlinear  case,  the  differential  equations  cannot  be  solved  for  the  sample  paths,  and 
frequency-domain  techniques  are  not  applicable. 

The  noise  sources  vi  and  V2  will  be  Thevenin  form  voltage  sources  with  power  spectral 
density  2kTR,  one  placed  in  series  with  each  resistor.  The  open-circuit  differential  equations 
become 

dq  q  1 
~dt  ^ 

and 


d(j) 

dt 


R/) 

~T~ 


+  VN2- 


(37) 


These  are  stochastic  differential  equations.  Rather  than  solving  the  differential  equations 
explicitly  for  sample  paths  and  then  taking  expectations,  the  expected  paths  can  be  solved 
for  directly.  Since  expectation  and  differentiation  are  linear  operations,  their  order  may  be 
interchanged,  resulting  in  differential  equations  for  the  means. 


and 

^  ^  {-f}  +  ^ 

These  simple  linear  differential  equations  can  be  solved  explicitly: 

E  {g(r)}  =  E  {g(t)}  e-h-b/(RC)^  ^ 

and 

E  {(/{t  +  t)}  =  E  {(/{t)}  T>t.  (39) 

The  variance  will  also  be  used  in  the  calculations.  The  quickest  way  to  obtain  the  variance 
is  to  appeal  to  the  equipartition  theorem,  which  says  to  each  degree  of  freedom  corresponds 
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an  energy  of  \kT.  The  capacitor  charge  is  a  degree  of  freedom,  so  the  capacitor  will  have 
an  expected  stored  energy  at  equilibrium 


E  {^]  =  l^T  ^  E  {fit}}  =  CkT, 


(40) 


and  the  expected  energy  stored  in  the  inductor  is 


^  E  =  LkT. 


(41) 


Because  the  noise  processes  and  stationary,  these  variances  have  hxed  values  for 

all  t. 

The  open-circuit  output  voltage  V  can  be  expressed  as 


.r  _  q  d4>  _  q  R4> 

dt  C  L 


(42) 


so  the  autocorrelation,  for  r  >  t,  is 

Rvvit-r)  =  E{V{t)V{T)} 


=  E 


q{t)  _ 

C  L 


EVN^it) 


q{r)  R(t){r) 


C 


i£{,(l),(r)}- As 

-^E{4,{t)q(r)]  +  ^^E 


L 

1 


+  njv2(r) 


0(^)}  +  ■^E{q{R)  ^'iV2(r)} 

0(^)}  -  VN^ir) 

+  g(^)}  -  0(r)|  +  E^^VN2{t)  njv2(r)| 

t)  VN^ir) 

jE^^VN^it)  0(r)|  +  E^^VN^it)  n7V2(r)|. 


R 


4>(t)  \  -  jE 


(43) 


The  last  equality  follows  because,  since  and  are  independent,  g(-)  is  uncorrelated 
with  either  (;/)(■)  or  nAr2(-)-  Further,  for  r  >  t,  because  of  the  independent  increments  property 
of  the  noise  process,  (pit)  is  uncorrelated  with  VN2iT),  so 


E  {(pit)  VN2iT)}  =  0,  T>t. 
For  the  Gaussian  white- noise  process  VN2it), 


-E'|n7V2(t)  'yAr2('r)|  =  2kTR6it  —  r). 


(44) 


(45) 
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Using  the  means  and  variances  found  above,  the  other  expectations  are  easily  computed  for 

T  >  t. 


=  E 
=  E 
=  E 


=  UfcT  T>t 


(46) 


E 


0(r)}  =  E|0(f)E|0(r)  0(t)}} 

{0(f)  [0(f)e-(^-*)^/^]} 


U|0(f)  [0(f)e-("-*)-f^/^l 

E|0(f)  0(f)} 

LkT  r>0 


(47) 


The  last  remaining  term  is  a  little  more  tricky  to  calculate.  Starting  with  the  differential 

(48) 


equation  (37)  and  multiplying  through  by  0(r)  yields 

d(f){t) 


0('r)  =  +  '^^2(^) 


dt  L 

Take  expectations  on  both  sides  of  the  equation, 

and  then  interchange  the  order  of  differentiation  and  expectation 

^(t)}  =  -^E  +  E  . 


(49) 


(50) 


The  term  E  {0(f)  0(t)}  was  just  calculated  in  Eq.  (47).  Using  it  and  its  derivative. 


yLkT  =  -y  {LkT  +  E  {vnM  0(r)}  ,  r  >  f. 


or  equivalently. 


E\vN,(t)  =  2RkT  t  >  i. 


(51) 


Now  that  the  expectations  have  been  calculated,  Eqs.  (45),  (46),  (47),  and  (51)  can  be 


substituted  into  Eq.  (43), 


Ryy{t  -  r)  =  ^  {CkT  ^  {LkT  e-("-b«A)  _  ^  {2RkT 


C2\  ^  '  L2 

+2kTR6{t  —  T),  T>t. 


(52) 
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Under  the  matching  condition  x  ~  time  constants  are  equal,  and  further,  the 

coefficients  are  related  by 

kT  _  kTR^ 

~C  ~  L  ■ 

Thus,  the  autocorrelation  simplihes  to 

Rvv{i  —  t)  = ‘^kTR  —  t),  T>t.  (53) 


For  T  <  t,  the  calculations  must  be  redone.  In  this  linear  case,  it  is  simple  to  re-express 
the  differential  equations  for  the  means,  under  the  condition  t  <  t.  For  example,  Eq.  (38) 
becomes 

E{g(f)}  =  T<t.  (54) 

Proceeding  along  these  lines  will  show  that  the  expressions  (46 )  and  should  have  —  |r  — 1| 
rather  than  —  (r  —  f)  in  the  exponents.  Such  is  not  the  case  for  (|51[),  which  vanishes  for 
r  <  f  by  the  independent  increments  property.  However,  the  corresponding  term  is  found  in 


E  vn2{t)},  which  no  longer  vanishes  as  asserted  in  (44),  but  instead  yields 


E 


VNo 


(r))  =  T<t. 


(55) 


In  the  nonlinear  case,  a  more  complicated  reversibility  argument  will  be  necessary  to  calculate 
this  last  expectation. 
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6  Nonlinear,  Noise-Free  Case 


In  this  section,  the  capacitor  and  inductor  are  allowed  to  be  nonlinear,  as  shown  in  Figure 
The  functions  /(■)  and  h(-)  are  assumed  to  satisfy  the  following  standard  assumptions  [8l[9l 

[ID]: 

.  /(O)  =  0,  h{0)  =  0 

•  h  and  /  are  continuously  differentiable  functions,  and  for  all  values  of  the  arguments 
and  some  fixed  e  >  0, 


—  >  e  >  0  and 
dq 


dh 

d(j) 


>  e  >  0. 


The  function  h{-)  need  not  be  odd.  This  system  is  a  special  case  of  the  situation  considered 
in  Ref.  [10],  where  the  inductors  were  required  to  have  odd  characteristics  in  order  that  the 
flux  random  processes  be  reversible.  The  system  matrix  is  diagonal  and  hence  symmetric 
because  the  states  are  decoupled  for  any  driving  voltage  V  (that  does  not  depend  on  I). 

The  circuit  differential  equations  are 


(k  ^  V  -  /(g) 
dt  R 

ft  = 

and  the  port  current  is 

The  goal  again  is  to  show  that  I  =  V/R,  which  requires 


(56) 

(57) 

(58) 


fig) 

R  ' 


(59) 
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This  is  not  a  sufficient  condition,  because  as  yet,  there  is  no  relation  between  q  and  (j). 

The  two  parallel  branches  can  be  used  to  write  two  distinct  expressions  for  the  port 
voltage  V: 


=  nh{4>)  +  ^ 

=  «§  +  /(«)■ 

If  fi^)  =  h{(j)),  then  these  last  two  expressions  combine  to  yield 

dq  1  d(j) 
dt  R  dt 

With  zero  initial  conditions  for  q  and  0,  it  follows  that 


(60) 

(61) 


(62) 


9(1)  = 


and  therefore  the  full  matching  condition  is 


4>(t) 

R  ’ 


h(0) 


or  /(g)  =  R  h{Rq). 


(63) 


(64) 


This  condition  is  really  a  constraint  on  the  functions  /  and  h,  and  not  on  the  dummy 
variables  g  or  It  may  be  more  helpful  to  express  it  as 


f{x)  =  R  h{Rx). 

Alternately,  if  one  considers  f'{q)  =  1/C(g)  as  the  reciprocal  of  the  incremental  capacitance 
and  h//)  =  1/L((/)  as  the  reciprocal  of  the  incremental  inductance,  then 

L(0)  =  i?2c(g)|,=^/^,  (65) 

a  local  version  of  the  linear  matching  condition  L  =  R?C.  Under  the  full  matching  condition 
with  zero  initial  conditions,  the  deterministic  terminal  behavior  of  the  matched  nonlinear 
bridge  is  I  =  V/R  :  the  second-order  nonlinear  circuit  reduces  to  a  simple  linear  resistor. 
What  happens  under  non- zero  initial  conditions  for  g  and  /?  Observe  from  Eqs.  (56) 


and  (57)  that,  under  the  matching  condition  (64)  but  irrespective  of  U(-), 

d{(f)  —  Rq) 


dt 


=  -Rhicj>)+fiq) 

=  —  R  hi^cf))  +  R  h{Rq) 
=  -Rh{0{<P-Rq), 


(66) 


where  /  lies  between  /  and  Rq,  by  application  of  the  Mean  Value  Theorem.  It  follows  that 

|[0_i?g]2  =  -2Rh{0{(j)-Rqf 

<  -2Re  {(()- Rqf  ,  (67) 

using  the  assumed  lower  bound  on  the  derivative  of  h.  Therefore,  the  quantity  (0  —  Rq) 
decays  to  zero  at  least  exponentially  quickly.  After  a  sufficiently  long  time,  the  matched 
nonlinear  bridge  circuit  is  again  indistinguishable  from  a  linear  resistor. 
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7  Nonlinear  Noisy  Case 


For  this  section,  a  Norton-form  Nyquist-Johnson  noise  model  is  associated  with  each  resistor 
in  the  circnit,  as  in  Figure  We  would  like  to  show  that  the  terminal  current  noise  of  the 
matched  bridge  is  the  same  as  that  for  a  single  linear  resistor.  To  first  order,  this  result  is 
clear.  Recall  that  the  incremental  capacitance  and  inductance  satisfy  T(0)  =  C{q)\q=fj,iR. 
A  linearization  about  the  noise-free  equilibrium  operating  point  (g,  0)  for  a  d.c.  applied 
voltage  of  a  nonlinear  matched  circuit  will  yield  a  matched  linear  circuit.  By  superposition, 
the  noise  current  for  the  linearized  circuit  is  unaffected  by  the  applied  voltage.  The  point 
of  this  section  is  to  show  that  this  equivalence  holds  exactly,  even  for  high  temperatures  or 
strong  nonlinearities  for  which  the  noise  could  drive  the  circuit  out  of  the  valid  region  of 
linearization. 

The  circuit  is  described  by  stochastic  differential  equations  (SDE’s): 


dq  _  V  -  f{q)  . 
dt  R 


(68) 


d(j) 

dt 


V  —  R  h{(f))  —  Rijsi,^, 


(69) 


where  and  1^2  are  independent  Gaussian  white  noise  processes  with  power  spectral  density 
2kT/R.  The  port  current  is 


I  = 


, , ,,  ^  r  /(,) 


tN2 


(70) 


One  might  be  tempted  to  use  the  matching  condition  (59)  and  immediately  conclude  I  = 


V/R  —  iN^.  However,  the  full  matching  condition  (64)  does  not  hold,  because  the  derivation 


was  for  a  different  excitation:  q  and  0  no  longer  satisfy  q  =  (j)/R,  because  they  are  now 
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driven  by  independent  noise  sources.  So,  the  power  spectrum  of  /  must  be  calculated  in  a 
more  methodical  way. 

Before  proceeding  to  study  the  noise  power  spectrum,  we  show  that  the  nonlinear  inductor 
and  capacitor  cannot  “rectify”  the  noise.  Rectification  would  cause  incorrect  “average” 
behavior,  or  first-order  statistics  of  the  circuit,  such  that  it  would  be  pointless  to  study  the 
second-order  statistic  of  the  power  spectral  density.  We  will  show  that 


E{I}  =  V/R 

for  any  (deterministic)  V;  then  specihcally  for  R  =  0,  /  will  be  zero-mean. 
Taking  expectations  on  both  sides  of  Eq.  (70), 

B{/}  =  B{fcW}  +  -  Whh  _  0. 


(71) 


(72) 


The  expectations  of  f{q)  and  h{(f))  may  be  computed  by  integration  against  the  densities  for 
q  and  0.  For  a  random  variable  driven  by  Gaussian  white  noise,  the  Fokker-Planck  equation 
(FPE)  [m  [12]  converts  the  stochastic  differential  equation  into  a  partial  differential  equation 
in  the  density  of  the  random  variable.  For  the  circuit  of  this  section,  the  densities  obey 


dq 

dt 

dt 


V-f{q) 

R 

Rh{(j)) 


—  iNo 


R 


dpq 

dt 

dt 


dq 

A 

dcf) 


V-f{q) 

R 


-Pq 


+ 


kT  d^Pq 


R  dq^ 


V-Rh{ct>))p^ 


+  kTR 


d(j)^ 


(73) 

(74) 


Using  the  matching  condition  /(g)  =  R  h{Rq),  these  two  equations  become  identical  up  to 
a  scaling.  Suppose  a  density  p^{(l),t)  satisfies  Eq.  (74).  Then  use  the  change-of- variables 
formula  for  probability  densities  to  define 


pq{t,  q)  =  p<p  (U  0(g))  =  p<p{t,  Rq)  R 


and  substitute  this  dehnition  into  (73)  along  with  the  matching  condition  (59): 

kT  d^ 


V  -  Rh{Rq) 
R 


Rpd^it,  Rq) 


+ 


R  dq^  - 


Rp^{t,Rq) 


(75) 


(76) 


One  may  divide  by  R  to  remove  that  factor  from  the  left-hand  side,  then  group  a  factor  of 
R  into  each  d/dq\ 


Rq)  = 


d 


d{qR) 


V  -  Rh{Rq) 
R 


P<t>{t,Rq) 


d^ 


R  d{qR) 


;P^{t,Rq).  (77) 


The  variables  g  and  /  in  the  Fokker-Planck  equations  are  dummy  variables;  the  change  of 
variables  according  to  0  =  qR  does  not  assert  any  relationship  between  the  waveforms  q{t) 
and  0(f)  in  the  actual  circuit.  We  have  just  transformed  Eq.  (73)  into  Eq.  (74).  The  densities 
corresponding  to  zero  initial  conditions  (delta  functions)  also  satisfy  Eq.  (75)  at  f  =  0: 


5(g)  =  RS{cl))\^=Rq 


(78) 
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The  densities 


Pq(9) 

pM 


Aq  exp 

1  1 

^\v-f{q))dq 

0 

■  1 

fv 

- 

A^  exp 

j 

(70 

(79) 

(80) 


where  Aq  and  A(j)  are  normalization  constants,  are  the  steady-state  solntions  to  Eqs.  (73)  and 
(74),  and  also  satisfy  Eq.  (75)  at  t  =  0.  (See  Appendix  II  for  a  derivation  of  these  densities 


and  a  discnssion  on  the  meaning  of  steady-state.) 
satisfy  Eq.  (75)  for  all  time.  Therefore, 


Thns,  the  solntions  of  Eqs.  ( 73 )  and  ( 74 ) 


E{fiQ)}  = 


r*H-co  /*+oo 

f{x)pq{x)  dx  =  /  R  h{Rx)  p^{R  x)  R  dx 


r’  +  CO 


=  R 


Kv)  p<i>{y)  dy 


=  RE{hm, 


(81) 


where  the  dnmmy  variable  y  =  Rx  (with  dy  =  R  dx)  was  nsed  to  go  between  the  hrst  and 
second  lines.  Therefore,  the  relationship  E{I}  =  V/R  has  been  established. 

The  short-circnit  {V  =  0)  power  spectral  density  is  the  Fonrier  transform  of  the  antocor- 
relation  of  /,  denoted  Ru-  For  V  =  0  and  r  >  t, 


Riiit  -t)  =  E 


^(0(^))  -  -  iN2{t)]  \h{(j){T))  - 


R 


R 


E{hm))  mr))}RE{l^  ^ 


(82) 


The  perhaps  snrprising  simplihcation  nnder  which  Eve  of  the  nine  terms  generated  by  the 
prodnct  disappear  occnrs  for  the  same  two  reasons  as  in  the  linear  case.  Firstly,  for  17  =  0, 
the  two  circnit  branches  are  disjoint,  so  that  h(0(-))  is  not  correlated  with  f{q{-))  or 
regardless  of  the  time  argnments.  Secondly,  f{qit))  depends  only  on  past  valnes  of  iN2(t) 
and  hence  is  nncorrelated  with  *^2(1')  for  r  >  thy  cansality  and  the  independent-increments 
property  of  iv2(‘)- 

Using  the  relation  between  the  densities  for  q  and  0  again. 


E  {f{q{t))  /(g(r))}  =  R^E  {h(0(f))  h(0(r))}  ,  (83) 


even  thongh  the  expectations  have  not  been  calcnlated  explicitly  (see  Appendix  III,  however, 
for  an  explicit  verihcation  that  this  identity  holds).  Therefore,  if  it  can  be  shown  that 

E{f{q{t))  f{q{T))}  =  E{iN2{t)  /(g(r))},  (84) 

this  wonld  leave  E  {iN2it)  *W2('r)},  which  is  exactly  the  desired  antocorrelation. 
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(85) 


In  the  V  =  0  situation,  the  stochastic  differential  equation  (SDE)  for  q  becomes 

dgi.t)  /(lit))  ,  , , 

If  we  multiply  through  by  /(q'(r))  and  take  expectations, 

^  f{q{r))y  (86) 

The  dummy  time  indices  t  and  r  may  be  interchanged,  corresponding  to  writing  the  SDE  in 
T  and  multiplying  through  by  f{qit)),  to  get 

/(?(^))|  =  -^El^fiqiT))  /(g(t))}  -  E|i^2(r)  /(g(t))}.  (87) 

By  causality  and  independent-increments,  for  t  >  t, 

E{iN2{r)  f{qit))}  =  0. 

Since  differentiation  and  expectation  are  linear  operations,  their  order  can  be  interchanged. 

^  ^  /(?(^))}  =  /(9(^))} 

Further,  since  the  stochastic  process  is  stationary,  the  process  q(t)  is  also,  and  thus  the 
expectation  is  only  a  function  of  the  difference  {t  —  t).  Let  F(-)  be  defined  by 


F(t-r)  =  E{q{t)  /(g(r))} 


E{q{T)  f{q{t))}  =  F{T-t). 


As  a  consequence  of  the  assumptions  on  /(■)  and  Eq.  (85),  q{t)  is  a  reversible  process 
US].  This  means  that  for  all  ti  and  t2, 

Pr  [a  <  q{ti)  <  a  +  da,  f3  <  q{t2)  <  /9  +  d/9]  ='Pi[j3  <  q{ti)  <13  +  d/9,  a  <  q{t2)  <  a  +  da] 
As  a  consequence  of  this  reversibility,  F  is  an  even  function: 

F(t-r)  =  E{q{t)  f{q{T))}  =  [ [  af{b)p{q{t)  =  a,q{T)  =  b)  dadb 


(88) 


af{b)  p{q{T)  =  a,  q{t)  =  b)  da  db  =  E  {g(r)  f{q{t))} 
=  F{T-t), 


where  p{- ,  ■)  represents  the  joint  probability  density  of  its  two  arguments,  and  equality 
between  the  first  and  second  lines  follows  from  reversibility.  Since  F{-)  is  an  even  function, 
F'(-)  must  be  odd.  Hence,  F\t  —  t)  =  —F'{t  —  t),  and 


d 


d 


■:nE{q{t)  f{q{T))}  =  —F{t-T)  =  F\t-T) 


dt 


dt 


d 


-F(T-t)  =  --F(T-t) 

-<E{q(T)S(qm- 


(89) 
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Thus,  by  adding  together  (86)  and  (87),  the  left-hand  sides  cancel,  so  that 


0  =  -^E{f{q{t))  /(g(r))}  -  ^E{f{q{T))  f{q{t))}  -  E  {iN^{t)  /(g(r))}. 


or,  equivalently. 


^  E  {/(9(*))/(9(t))}  =  -E  {!«,(()  /(?(r))}  , 


(90) 


which  differs  from  (84)  by  only  a  common  factor  of  R.  Therefore, 


2kT 

Riiit-r)  =  E{I{t)  I{t)}  =  E  “  ^)- 


The  nonlinear  noisy  matched  bridge  circuit  has  a  short-circuit  noise  current  spectral  density 
precisely  the  same  as  a  single  linear  resistor.  Appendix  IV  will  show  that  the  noise  process 
is  Gaussian  white  noise. 
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8  Nonequilibrium  Situations 

The  noise  results  so  far  have  all  been  equilibrium  results.  For  nonequilibrium  situations,  the 
noise  component  of  the  total  current  is 

n{t)  =  I{t)- 

As  mentioned  in  the  introduction,  there  are  some  situations  in  which  the  noise  current  of 
the  matched  bridge  is  not  statistically  equivalent  to  the  noise  of  a  single  linear  resistor. 
However,  the  result  does  hold  in  some  specihc  nonequilibrium  situations.  For  example,  if  the 
two  resistors  of  the  matched  linear  bridge  are  maintained  at  the  same  constant  temperature 
by  a  large  thermal  reservoir,  the  noise  component  of  the  output  current  will  be  Gaussian 
white  noise  with  power  spectral  density  2kT / R,  for  any  driving  voltage  V{t).  This  is  a 
result  of  superposition,  in  which  the  deterministic  driven  response  is  decoupled  from  the 
noise  response.  Since  the  proper  behavior  has  been  proven  for  the  responses  individually, 
the  total  response  will  be  correct.  The  explicit  calculations  are  found  in  Appendix  I. 

In  the  nonlinear  case,  we  cannot  use  superposition  to  establish  our  result.  Nevertheless, 
it  is  interesting  that  the  result  will  still  hold  if  the  driving  voltage  is  constant.  Unfortunately, 
the  result  does  not  hold  in  a  time- varying  circuit,  and  this  failure  casts  doubts  on  the  hopes 
of  establishing  the  general  nonlinear  nonequilibrium  result  for  time-varying  voltage. 


8.1  Non-zero  drive  voltage 

Consider  the  circuit  of  Figure]^  again  described  by  Eqs.  (68)  to  (70)  with  the  matching 
condition  (59)  in  force,  but  now  with  an  arbitrary  driving  voltage  applied  to  its  terminals. 
Denote  by  Rnn  the  autocorrelation  of  n{t)  =  I (t)  —  V{t)/R.  Then 


9kT 

Rnnit-T)  =  —6it-T), 


(91) 


provided  that  either  (a)  the  circuit  is  in  fact  linear,  or  (b)  the  voltage  V{t)  is  constant. 
For  an  arbitrary  driving  voltage  V{t),  the  autocorrelation  of  n{t)  is 

Rnn{t,T)  =  E  I  [h(0(t))  -  [h(0(r))  - 


=  E 


R 

h{<P{T))]-E 


/(g(^)) 

R 


R 

-E 


-e\1^  h(0(x))|  +  e{1^  l^]  +  e\1^ 


I  R 

-EliN^it)  h(0(r))  \  +E\iN^{t) 


I  R 

+  E\iN^{t)  iN^ir) 


R 

As  in  the  short-circuit  case,  h(0(-))  is  not  correlated  with  iN2{'),  so  that 


;{h(0(t))  iiV2(^)}  =  ^{*7V2(^)}  =  0. 
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The  process  h(0(-))  is  also  uncorrelated  with  f{q{-)),  but  these  processes  are  no  longer 
zero-mean,  so  we  can  only  separate  the  expectation  of  the  product  into  the  product  of  the 
expectations, 


Because  of  the  similarity  of  the  Fokker-Planck  equations  (73)  and  (74), 


E 


1 


for  all  times  t  (or  r),  and  also 

E 

The  autocorrelation  can  thus  be  simplihed  to 


h(4,(r))]  =  £{h|i)  fjm 


R 


Rnnit^r)  =  2E 
+E 


R 


R 


-2E 


R 


*iV2(^)  \  +E<iN2{t) 


R  ) 


R 


R 


+  E\iM^{t)  (92) 


Using  the  trick  of  Section  7,  multiplying  both  sides  of  the  differential  equation  for  q{t)  by 
/(g(r)),  we  obtain 


I- 


{f{(l{r))  g(t)}  =  s|/(g(r))|  -  E|/(g(t))  /(g(r))|  -  El^iN,{t)  /(g(r))|.  (93) 


Define 


F(f,r)  =  E|/(g(r))  g(f)|, 
SO  that  the  autocorrelation  may  be  expressed 


TfT--® 


1  — Jin  — 


//(9W) 


V(t) 


R 


-E 


/(g(r)) 

R 


E 


dE(t,  t)  ^  dE(t,  t) 


dt 


dr 


+  Eli  N^{t)  In  At) 


mt))\ 

(94) 


For  arbitrary  time- varying  V{t),  no  further  simplihcation  is  apparent.  However,  if  V  is 
constant,  he.,  the  system  is  at  steady-state,  then 

WiP"' 


and  by  taking  expectations  of  both  sides  of  the  differential  equation  (68), 

(V{t)-f{q{t))\ 


R 


E{iN2{t)}^ 
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so  that 


V{t)  =  E< 

Further,  q{t)  is  again  a  stationary  random  process,  so  that  F{t,T)  =  F{t  —  r)  and 


dF{t,  t)  ^  dF{t,  t) 


dt 


dr 


and  therefore 


Rnn{t,T)  =  E<  (t)  \  =  ——  -  t). 


=  0, 


2kT 


R 


8.2  Linear  time- varying  elements 


Up  to  this  point,  every  case  that  we  have  been  able  to  solve  has  yielded  the  correct  result,  and 
we  cannot  conclude  anything  from  our  failure  to  simplify  the  equations  for  the  nonlinear  case 
with  time-varying  drive  voltage.  The  following  case  will  hnd  our  hrst  nontrivial  exception. 

Suppose  the  energy  storage  elements  in  Figure  were  linear,  but  time  varying.  It  is 
sufficient  to  consider  the  short-circuit  (undriven)  behavior.  The  circuit  differential  equations 
are 


(95) 

(96) 


and  the  port  current  is 


dcj) 

R^it) 

dt 

Lit) 

dq 

(lit) 

dt 

RCit) 

T(t\  — 

(j){t) 

(lit) 

1  {i) 

Lit) 

RCit) 

condition  is 

of  course 

iN2{R)- 


L{t)  =  C{t). 


(97) 


(98) 


The  differential  equation  for  q{t)  can  be  solved  explicitly  in  terms  of  sample  paths  of  the 
noise  process  iN2(t): 


q(t)  =  exp 


ds 


g(0) 


lN2(a)  exp 


ds 


Jo  RC(s) 

The  autocorrelation  function  for  the  port  current  for  r  >  f  is 


RC{s) 


da 


(99) 


Riiit,T)  =  E 


q{t) 


Lit)  RCit) 
1 


i?2C(t)C'(r) 
+E\it^^it)  iisi^ir 


-  iN2  {t) 

+ 


Xr)  qir) 


_L(r)  RCir) 
1 


-  ( r ) 


L{t)L{T 


-E 


0(r)  \  + 


RCir) 
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where  the  other  terms  vanish  because  the  variables  are  uncorrelated  (and  zero-mean)  or  by 
causality,  as  argued  previously.  Again  by  appeal  to  the  Fokker-Planck  equations  and  the 
matching  condition  (98),  it  can  be  shown  that 

1 


C{t)  C{t) 


Q{r)  = 


m  Lir] 


E 


0(r)|. 


Thus,  in  order  that  the  short-circuit  current  noise  have  the  proper  autocorrelation,  it  must 
be  shown  that 


C{t)  C{t) 


(100) 


Two  quick  calculations  from  Eq.  (99)  yield 

/■*  ds 


E{qit)qiT)}  =  exp 


'0 


RC{s) 


exp 


RC{s) 


X  E{g2(0)}  + 


2kT 


exp 


ds 


ds 


and 


2kT 

E{q{T)  iN2{t)}  = - 


ds 


RC{s) 


Substituting  these  into  Eq.  (100),  the  test  is 

/■*  ds 


7 

0  = 


exp 


R^C{t)C{T 

x(E{q\0)}  + 


RC{s) 
2kT 


exp 

rt 


2kT 


RC{t)  R 
/■*  1 


exp 


R 

ds 


exp 


'0 


RC{s] 

ds 


exp 


■ds 


RC{s) 


ds 


da 


RC{s) 


(101) 

(102) 


RC{s) 


exp 


RC{s) 
ds 


>0 


da 


=  exp 


'0 


RC{s) 


ds 


E{q\0)}  + 


2kT 


RC{s) 

/  exp 

'o 


ds 


RC{s)_ 


da 


RC{t)  2kT 
2  R 


exp 


ds 


L-'o 


RC{s) 


=  E{q\0)}  + 


2kT 


exp 


hr  1 

,  RC{t)  2kT 

[Jo  RC{s)\ 

aa  2  R 

L“7o  RC{s)\ 

(.103) 


Differentiating  by  t  will  yield  a  necessary  condition  for  the  equation  to  be  true: 


V  2kT 

0  =  -^exp 


ds 


RC{s) 


,  dC(t)R2kT 
da - t: — ^  exp 


RC{t)  2kT 
2 

dC{t) 


exp 


ds 


dt 


kT  exp 


'0 


'o  RC{s) 
'■*  ds 
RC{s) 


dt  2  R 

imt)) 


ds 


RC{s) 


26 


Thus,  C  must  be  a  constant.  In  this  case,  the  integrals  in  Eq.  (103)  can  be  computed,  and 
if  the  system  starts  at  equilibrium,  he.,  E{q‘^{0)}  =  kTC,  then  this  condition  is  sufficient 
as  well  as  necessary.  Of  course,  if  O  is  a  constant,  then  the  bridge  is  simply  the  linear, 
time-invariant  case  that  was  already  considered  in  Section  3.  A  time-varying  bridge  circuit 
can  be  distinguished  from  a  single  linear  resistor,  although  it  is  hard  to  imagine  the  black 
box  of  the  layman’s  test  corresponding  specifically  to  this  case. 

The  importance  of  this  case  comes  from  the  following  analysis:  For  a  driving  voltage 
V{t)  signihcantly  larger  than  the  noise,  one  could  solve  the  deterministic  system  and  then 
compute  an  approximation  for  the  noise  behavior  by  linearization  about  this  time-varying 
solution.  This  approximation  would  behave  like  the  time-varying  linear  system  described 
above.  Since  the  second-order  statistics  for  that  system  are  incorrect,  we  believe  that  the 
second-order  statistics  for  the  nonlinear  system  driven  by  a  time-varying  voltage  will  not 
match  the  statistics  of  a  single  linear  resistor  driven  by  that  same  voltage. 
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9  Conclusions 


This  paper  has  considered  the  application  of  the  flnctnation-dissipation  theorem  (or,  in 
circnit  theory  terms,  a  result  relating  impedances  to  noise  spectra)  to  many  variations  of  a 
particular  bridge  circuit. 

The  linear  results  were  previously  known,  but  we  are  unaware  of  a  reference  that  in¬ 
cludes  all  of  the  special  cases  (the  circuit  and  its  dual)  and  different  approaches  to  the 
solution  (sample-path,  frequency-domain,  and  Fokker-Planck  equation  as  a  special  case  of 
the  nonlinear  results). 

Further,  we  have  indicated  an  extension  of  this  fluctuation-dissipation  theorem  to  a  non¬ 
linear  situation.  The  spectral  calculations  for  this  case  have  been  nontrivial,  calling  on  a 
reversibility  idea  and  martingale  theory. 

The  positive  results  hold  for  a  specihc  time-invariant  bridge  circuit,  linear  or  nonlinear, 
in  thermal  equilibrium  or  at  d.c.  steady-state.  The  linear  results  have  also  been  shown  to 
hold  for  the  dual  circuit  and  for  an  arbitrary  time-varying  drive  voltage. 

The  negative  results  in  Section  show  that  our  original  fluctuation-dissipation  conjecture 
is  not  correct  as  stated  and  must  be  limited  to  exclude  time- varying  networks  and  nonlinear 
networks  with  time-varying  inputs.  Is  the  modified  form  below  correct?  This  remains  an 
open  question  in  the  field,  and  some  of  the  ideas  in  [10]  may  be  of  assistance. 

Modified  Fluctuation- Dissipation  Conjecture  for  Circuits:  No  two  zero-state  de¬ 
terministically  equivalent  time-invariant  networks  can  be  distinguished  by  the 
terminal  noise  currents  at  any  d.c.  voltage  input  when  the  networks  are  in  sta¬ 
tistical  steady-state. 

The  assumptions  here  remain  those  in  the  paragraph  preceding  the  initial  formulation 
(see  the  Introduction),  including  LTI  Nyquist- Johnson  resistors  and  nonlinear  inductors  and 
capacitors.  Additional  assumptions  may  be  required  to  guarantee  reversibility  of  the  charge 
or  flux  random  processes. 

The  further  extensions  to  include  nonlinear  resistor  noise  models  or  multiterminal  circuits 
remain  completely  unexplored,  so  far  as  we  know.  (For  nonlinear  resistors,  progress  is 
hampered  by  lack  of  a  universally-accepted  noise  model  for  nonlinear  devices.) 


Appendix  I:  Linear  Noisy  Case  with  Non-zero  Drive 

The  result  of  Section  may  also  be  obtained  without  short-circuiting  the  terminals.  If  the 
terminals  of  the  linear  bridge  circuit  are  driven  by  a  deterministic  voltage  waveform  V{t), 
the  autocorrelation  of  the  output  current  must  reduce  to  the  autocorrelation  of  the  current 
driven  by  V (t)  applied  across  a  noisy  resistor,  namely, 


Rn{T)  =  E  {/(«)  I{t  +  r)}  =  hi*  ^  Hr)- 


R 


(104) 


This  result  can  be  claimed  immediately  as  a  result  of  superposition,  since  in  this  linear  case, 
the  deterministic  response  and  the  noise  behavior  are  decoupled.  Some  calculations  verify 
that  claim. 

The  solutions  to  the  differential  equations  ([^  and  ([^  with  drive  terms  are 


Qit)  = 

.it)  =  j 

The  autocorrelation  is 
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Rii{r)  =  E 
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+  E{iN2it)  iN2it  +  T)}. 
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Working  backwards,  it  will  be  shown  that  all  terms  bnt  the  two  in  the  last  line  vanish  or 
cancel.  Of  conrse,  E  {iN2(t)}  =  0,  removing  two  terms.  A  qnick  calcnlation  nsing  Eq.  (106) 
shows  that  i(-)  is  nncorrelated  with  even  for  E  7^  0. 

The  single  term  expectations  may  be  computed  from  the  pathwise  solution  ( 105 )  for  q: 

rt 
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and  (106)  for  i: 
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Therefore,  under  the  matching  condition  R/L  =  1/RC, 
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(109) 


and  the  same  holds  if  the  time  arguments  are  switched.  (Note  that  E{i{t)}  7^  E{i{t  +  r)} 
in  general,  but  the  cancellation  does  not  require  this.)  This  leaves  a  more  manageable  set  of 
terms  in  the  autocorrelation. 

Rjj{t)  =  E{i(t)i(f  +  r)}-^E{z(f)g(f  +  r)}-^E{g(t)i(f  +  r)} 

+  E  {q{t)  q{t  +  t)}  +  ^  E  {q{t)  +  t)}  +  ^  E  m (t)  q{t  +  t)} 


V{t)  V{t  +  t) 


+  E{iN2it)  iN2it  +  T)}. 


(110) 


Then,  by  using  the  pathwise  solutions  (105)  and  (106),  the  remaining  expectations  can  be 
calculated  just  as  they  were  in  the  previous  section. 
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Under  the  matching  condition,  ^  C,  and  R/L  =  l/i?U,  so  that 


E  {i{t)  i{t  +  r)}  =  {q{t)  q{t  +  r)}  . 


The  expectation 


E{iN2{t)  q{t  +  T)}  = 
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,  r  >  0, 
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(114) 


found  in  Eq.  (01  is  still  valid  because  U(-)  is  uncorrelated  with  iN2{')-  However,  the  two 
remaining  terms  of  Eq.  (110)  were  not  calculated  in  the  previous  case  because  g(-)  and  i(-) 
were  uncorrelated.  This  is  no  longer  true  with  a  non-zero  applied  voltage;  instead. 
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because  the  deterministic  voltage  driving  each  process  is  the  same  (that  is,  completely  corre¬ 
lated),  even  though  the  random  current  processes  driving  them  are  independent.  Since  the 
time  constants  are  equal  under  the  matching  conditions. 
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(116) 
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The  term 


r-t+r 


^-{t-t')R/L  ^-{t+T-t")/RC  Y(f'^  V{t'‘ 


'  — oo  J  —oo 


appears  in  four  terms  of  Eq.  110:  twice  positively  {E  {i(t)  i(t  +  r)}  and  E  {q{t)  q{t  +  r)})  and 
twice  negatively  {E  {i(t)  q(t  +  r)}  and  E  {q(t)  i(t  +  t)}).  The  matching  condition  shows  that 
the  coefficients  are  equal  in  magnitude,  so  that  these  terms  drop  out  of  the  hnal  expression 
for  Rii{t). 

Similarly,  just  as  in  the  undriven  case,  the  term 


^  ^-\r\R/L  ^ 


i?2C'2 


kTC 


appears  in  four  terms:  positively  in  E  {i(t)  i(t  +  r)}  and  E  {q(t)  q{t  +  r)},  and  negatively 
but  doubled  in  E  q{t  +  t)}  for  r  >  0  or  E  +  t)  q{t)}  for  r  <  0.  Therefore, 

these  terms  again  cancel,  proving  the  claim  that 


Rii{r)  = 


V{t)V{t  +  T)  2kT 


i?2 


R 


5{t) 


(117) 
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Appendix  II:  Steady-State  Densities 


In  this  appendix,  we  are  searching  for  the  steady-state  densities  solving  the  Fokker-Planck 
equations  (73)  and  (74).  But  hrst,  a  word  about  what  steady-state  means,  and  what  the 
solutions  to  our  equations  look  like. 

In  the  absence  of  noise,  the  circuit  equations  for  the  nonlinear  bridge  circuit  are 


dq 

dt 

d(f) 

dt 


V  -  f{q) 

R 

V  -Rh{(j)). 


(118) 

(119) 


A  d.c.  equilibrium  or  equilibrium  operating  point  for  a  hxed  (d.c.)  voltage  V  is  any  point 
(g,  (j))  for  which  ^  ^  =  0.  Solutions  to  the  circuit  differential  equations  are  trajectories 

in  the  (g,  0)  plane. 

In  the  presence  of  noise,  the  circuit  equations  for  the  nonlinear  bridge  circuit  are  stochastic 
differential  equations: 


dq  _  V  -  f{q)  _  . 
dt  R 

^  =  V-Rh{(j))-RtN,. 
dt 


(120) 

(121) 


The  solutions  to  these  equations  can  be  expressed  in  two  different  ways.  Sample-path  so¬ 
lutions  are  trajectories  in  the  (g,  0)  plane  that  show  the  behavior  of  the  circuit  to  specific 
realizations  of  the  noise  processes  i^-^  and  ijy^.  Sometimes,  as  in  Section]^  the  trajectories 
may  be  computed  explicitly.  For  most  nonlinear  stochastic  differential  equations,  however, 
explicit  solutions  cannot  be  calculated. 

In  such  cases,  one  turns  instead  to  Fokker-Planck  equations  In  fact,  the  FPE  is 

also  frequently  used  in  the  linear  case,  when  one  is  interested  in  the  statistical  behavior  of 
the  circuit  rather  than  individual  sample  paths.  The  FPE’s  corresponding  to  the  stochastic 
differential  equations  above  are,  respectively. 


dpg 

d 

V-fiq) 

dt 

dq 

R 

dp^ 

dt 

d 

d(j) 

(v  -Rh 

+ 


kT 

R  dq^ 


P(j> 


+  kTR 


d(j)^ 


(122) 

(123) 


The  solution  to  an  FPE  is  a  probability  density  p  which  lives  in  the  infinite-dimensional 
space  of  real-valued  nonnegative  smooth  functions  (a  density  is  also  normalized  such  that 
its  integral  over  all  space  is  1). 

A  steady-state  density  satishes  ^  =  0.  Thermal  equilibrium  (not  to  be  confused  with  a 
d.c.  equilibrium  point)  for  this  circuit  is  the  steady  state  with  V  =  0. 

We  are  now  prepared  to  search  for  the  steady-state  solutions  to  the  FPE’s  above  for  a 
constant  V.  Starting  with  the  hrst  equation,  for  the  probability  density  of  the  capacitor 
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charge, 


0 


dq 

dq 


V-f{q) 

R 

V  -  f{q) 

R 


+ 


-Pq 


kT  d^Pg 
R  dq^ 

R  dq  ) 


(124) 

(125) 


The  density  p  must  decay  to  zero  at  q  =  ±cx)  so  that  p  is  integrable;  its  derivative  dp/dq 
hence  necessarily  also  decays  to  zero.  Therefore,  at  steady-state,  not  only  must  the  inside 
of  the  large  parentheses  be  a  constant  with  respect  to  g,  that  constant  must  be  zero.  (For 
multidimensional  systems,  one  must  require  reversibility  or  detailed  balance,  which  are  deep 
physical  concepts  that  essentially  say  that  one  cannot  tell  if  time  is  running  forwards  or 
backwards.  One  could  tell  the  direction  of  time  if  probability  were  exiting  ai  q  =  oo  but 
entering  at  g  =  — cx)  to  maintain  a  total  probability  of  unity.)  Thus,  we  are  looking  only  for 
solutions  to 


0 


V-f{q) 

R 


-Pq 


kT^\ 
R  dq  )  ' 


or,  equivalently. 


dpq  _  V  -  /(g) 
dq  kT 


(126) 


The  solution  is 


Pq{q)  =  Agexp 


(127) 


where  Ag  is  a  normalization  constant. 

Similarly,  solving  for  the  probability  density  of  the  inductor  flux  from 


d(j)  L 


reduces  to 


V-Rh{^)^p^ 

d^  ^  V-R  h{(j)) 
d4> 


RkT 


+  kTR 


P(j>- 


d(jA 


The  steady-state  solution  is  then 

p<^(0)  =  A^exp 


(128) 


(129) 


(130) 


where  A^p  is  a  normalization  constants. 
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Appendix  III:  Explicit  Covariance  Calculation 

The  identity  to  be  proven  in  this  section  is 

E  /(g(r))}  =  R^  E  {h(0(t))  h((/)(r))}  .  (131) 

In  earlier  drafts  of  this  paper,  the  identity  was  simply  asserted  to  follow  from  the  Fokker- 
Planck  eqnations  for  the  densities  of  q  and  0,  Eqs.  (73)  and  (74),  respectively.  Fnrther,  the 
FPE  is  really  a  conditional  density  evolntion  eqnation,  which  assnmes  some  initial  density 
that  then  evolves  in  time.  So,  it  shonld  not  trouble  us  that  we  shall  be  splitting  the  left-hand 
side  into 

E{f{q{t))  /(g(r))}  =  E|/(g(r))  E{/(g(t))|g(r)}} 
for  t  >  T.  We  shall  use  the  notation 

Pq{t)ix) 

for  the  probably  density  of  the  random  variable  q{t)  evaluated  at  x,  and 

iy\^) 

for  the  probably  density  of  the  random  variable  q{t)  evaluated  at  x,  conditional  on  g(r)  =  x. 
Then, 


E{f{q{t))f{q{T))} 

(1) 

(2) 

(3) 

(4) 

(5) 


=  E 


{fiqir))  E{f{q{t))\q{T)}^ 


b  +  OO 


'  — oo 
/•  +  00 


^+oo 


E? 

R^ 


fix)pq{T){x, 

R  h{Rx)pg(^r){x 
h{Rx)R 
h{x')  P4>{r){x') 
h{x')  P^{t){x') 


f{y)pq{t)\q{T){y\x)dy 


dx 


H-oo 


Rh{Ry)pq^t)\q{T){y\x)dy 


dx 


^+oo 


'  — OO 

r+oo 


'  — oo 
r+oo 


P  +  OO 


HRy)Pq{t)\qir){y\x)dy 


— oo 

H-oo 


d{Ry^Pq{t)\q(T)  I  y 


.J  — oo 
/•+00 


y' 

h{y  )Pq(t)\q(T)  (  ^ 


^  1  ^ 


dx 

dx' 

dx' 


where: 

(1)  is  writing  the  expectations  in  terms  of  densities 


(2)  uses  the  matching  condition  (64)  f{x)  =  R  h{Rx) 


(3)  uses  the  change-of- variables  (75)  for  single-time  densities,  Pq{x)  =  R  p^{Rx) 

(4)  rescales  according  to  x'  =  Rx 

(5)  rescales  according  to  y'  =  Ry 
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We  now  write  out  the  explicit  integral  for  the  right-hand  side  of  our  identity. 
E{h(0(t))  h(0(r))}  =  E{h(0(r))  E{h(0(t))|0(r)}} 


“+C» 


Kx)  P0(r)(a^') 


^+oo 


h{y')pm\<f>{r)  iy'W)  dy' 


dx 


So,  all  that  remains  to  showing  our  identity  is  to  show  that 


Pq{t)\q{T) 


X‘ 


1  P4>(tM('r)  (y  1^  ) 


If  we  start  with  the  FPE  (73)  expressed  in  the  dummy  variable  y, 

kT  d^pg{s,y) 


dpq{s,y)  ^ 

ds  dy 


v-f{y) 

R 


Pqi.s,y) 


+ 


R  dy^ 


(132) 


(133) 


we  can  compute  the  conditional  density  as  follows: 


Pqit)\q{T) 


dpq{s,y'/R) 

ds 


ds  +  Pq{T)\q(T)  (  ^ 


d 


:  d{y'/R) 

+5 

\R  R 


V  -  f{y'/R) 

R 


Pq{s,y'/R) 


kT  d‘^pq{s,y' /R) 
R  d{y'/Rf 


^  +  kTE^(s,y' /E)  ds 


Similarly,  from  Eq.  (74), 


P0(t)i0(T)  {y  ) 


;  dy' 

+R6{y'  -x') 

” -T^V  -Eh(y')]p,(s,y' IE) 
+R5{y'  -x') 

r  dp^{s,  y') 


+  kTR 


dyl2 


s,  y'  /  R)  ds 


(134) 


ds 

A 

dy'  L 


ds  P(^(r)|<^(T)  (?/  1^  ) 


V  -  Rh{y'))p^{s,y') 


+  kTR  ^  (s,  y')  ds 


du 


a 


-x^) 

The  initial  conditions  are  scaled. 


(135) 


Pq{T)\q(T) 


R 


X 


y  X 


R  R 


-^1  =d  (~5~  R  I  =  ^d{y'  -x')  =  Rp^(^rM^r)  {y'W) 


and  the  unconditioned  density  appears  linearly  in  the  integrals,  so  that  the  scaling  will  be 
preserved  as  we  evolve  away  from  the  initial  conditions. 
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Appendix  IV:  Martingale  Theory 


A  technical  question  remains:  whether  the  short-circuit  current  noise  process  is  in  fact 
Gaussian  white  noise,  or  whether  it  just  has  the  same  power  spectral  density.  Let  x{t)  be 
dehned  by 

x{t)  =  f  I{t)  dr. 

J  a 

Then  x{t)  is  a  martingale  on  any  finite  interval  a  <  t  <  b:  denoting  by  J^s  the  history  of  the 
process  up  to  time  s, 


E  {x{t)  \  Ts}  =  x{s)  +  E  <  /  I (t)  dr  >  =  x{s)  +  /  E  dr  =  x{s), 


because  /(r)  is  zero-mean  (for  V  =  0).  The  sample  paths  of  x{t)  are  continuous,  since  the 
right-hand  side  of  ( 70 )  consists  of  all  continuous  functions  na.  The  variance  of  x{t)  is  finite: 


E{x‘^{t)}  =  E 


I{t)  dr 


I{t')  dr 


t  rt 


E  dr  dr' 


*  2kT  ,  2kT  ^ 

o[t  —  T  )dT  dr  =  ——  [t  —  a)  <  oo. 


R 


R 


And  further,  for  t  >  s, 


E  <  {x{t)  —  x{s)y 


JF,  !>  =  E  \^x‘^{t)  —  2x{t)x{s)  +  x‘^{s)  \  Eg} 

E  {x\t)  \Es}-2E  {x{t)  I  Eg}  +  x\s) 

,  zTi  I  r  ^  I{t)  dr  [  I{T')dT' 


=  x^{s)  +  E 


t  rt 


E  {/(r)J(r')}  dr  dr  = 


t  nt 


'  s  J  s 

2kT 


2kT 


—  2x{s)x{s)  -\-  x^(s) 
6{t  —  T')dT  dr' 


R 


{t-s). 


The  process  x{t)  satisfies  all  of  the  hypotheses  of  Theorem  11.9  of  Ref.  [13],  from  which 
we  conclude  that  it  is  a  Brownian  motion.  Therefore  the  short-circuit  noise  current  current 
process  is  Gaussian  white  noise. 
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